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Abstract 

Let (G,P) be a bar framework of n vertices in general position in R d , for d < n — 1, 
where G is a (d + l)-lateration graph. In this paper, we present a constructive proof that 
(G,P) admits a positive semidefinite stress matrix with rank (n — d — 1). We also prove a 
similar result for a sensor network, where the graph consists of m(> d + 1) anchors. 

1 Introduction 

Let V(G) and E{G) be, respectively, the vertex set and the edge set of a simple edge-weighted 
graph G, where each edge has a positive weight dij. The graph realization problem (GRP) 
is the problem of determining whether there exists a realization of G in Euclidean space M. d , 
for a given dimension d. A (matrix) realization P = [p\, . . . , p n ] of G in M. d is a mapping 
P : V(G) — >■ M. dxn such that, if vertices i and j of 67 are adjacent, then the Euclidean distance 
between points pi € M. d and pj G M. d is equal to the prescribed weight dij on the edge We 
always assume that the points pi, . . . ,p n amnely span M. d . In order words, P is a realization of 
G if and only if p% , . . . , p n amnely span M. d and 

\\Pi ~ Pj\\ = dij for each edge € E(G). 

Throughout this paper, |jx|| denotes the 2-norm of a vector x. Also, we use to denote the matrix 
of all zeros of the appropriate dimension. See, e.g., [TZl UH1 IS El [T21 [T31 HH ISH [331 SI [T^l • 

The GRP and its variants arise from applications in various areas, such as molecular con- 
formation, dimensionality reduction, Euclidean ball packing, and more recently, wireless sensor 
network localization EJ EH Hfll E21 [20] . 
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Let P be a given realization of graph G with n vertices in Mr. A realization P together with 
a graph G is often referred to as a bar framework (or framework), and is denoted by (G,P). If 
P is the only realization of G in M. d , up to a rigid motion (e. g., translation or rotation), then 
we say that the framework (G, P) is globally rigid. However, if P is the only realization of G, up 
to a rigid motion, in all dimensions, then we say that the framework (G, P) is universally rigid. 

For a given framework (G, P) in M. d , define the (d + 1) x n matrix A such that 



A 



P 



(1) 



where e is the vector of all l's in W 1 . Matrices P and A are also respectively called the position 
matrix and the extended position matrix of the framework (G, P). The notion of a stress matrix 
plays a critical role in the characterization of the universal, as well as the global, rigidity of 
frameworks. An n x n symmetric matrix S is called a stress matrix of framework (G, P) if and 
only if 

AS = 0, (2) 

and 

S ij = 0,V(i,j)?E(G), (3) 

where A is the extended position matrix of (G,P). Note that the highest possible rank of a 
stress matrix S is (n — d — 1), and the zero matrix is a trivial stress matrix. 

The following theorem characterizes the universal rigidity of generic frameworks in terms 
of stress matrices. A framework (G, P) is said to be generic, or in generic position, if the 
coordinates of p\, . . . ,p n are algebraically independent over the integers, i.e., if there does not 
exist a non-zero polynomial / with integer coefficients such that f(pi, ■ ■ ■ ,p n ) = 0. 

Theorem 1. Let (G, P) be a framework of n vertices in generic position in M. d , d < n — 1. Then 
(G, P) is universally rigid if and only if there exists a stress matrix S of (G, P) such that S is 
positive semidefinite ( PSD ) and the rank of S is n — d — 1 . 

The "if" part of this theorem was proved independently in |13| and [1] , while the "only if 
part was proved in [15] . 

One of the major research topics in rigidity is whether a result similar to Theorem Q] holds 
if the assumption of a framework in generic position is replaced by the weaker assumption of a 
framework in general position. We say that framework (G, P) is in general position in M. d if no 
{d + 1) points of pi, . . . ,p n are affinely dependent. For example, points are in general position 
in M? if no 3 of them are collinear. It then easily follows that if framework (G, P) in W 1 is in 
general position, then every (d + 1) x (d + 1) square sub-matrix of the extended position matrix 
A, defined in ([1]), has rank (d + 1). Note that whether or not n rational points are in general 
position can be checked in time polynomial in n for any fixed dimension d, while the generic 
position condition is uncheckable. The following theorem, proved in [6] recently, shows that the 
"if part of Theorem [JJ still holds true under the general position assumption. 

Theorem 2. Let (G,P) be a framework of n vertices in general position in M. d , d < n — 1. 
Then (G, P) is universally rigid if there exists a stress matrix S of (G, P) such that S is positive 
semidefinite and the rank of S equals n — d — 1 . 
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However, it remains an open question whether or not the converse of Theorem [2] holds true. 
In this paper, we settle this question in the affirmative for frameworks (G, P) in general position 
when G is a (d+ l)-latertation graph. A graph of n vertices is called a (d + l)-lateration graph 
if there is a permutation ir of the vertices, vr(l), 7r(2), . . . , vr(n), such that 

• the first {d + 1) vertices, vr(l), . . . , 7r(d + 1), form a clique, and 

• each remaining vertex 7v(j), for j = (d + 2), . . . n, is adjacent to {d + 1) vertices in the set 
Ml),vr(2),...,vr(i-1)}. 

Such frameworks were shown to be universally rigid in [2D] and [23| . where several classes of 
universally rigid frameworks in general position were identified. 

In particular, we present a constructive proof that a framework (G, P) of n vertices in general 
position in M. d , where G is a (d + l)-lateration graph, admits a PSD and rank (n — d — 1) stress 
matrix S. We show that such a stress matrix S can be computed in strongly polynomial time, 
if the (d + l)-lateration ordering is known. We also show that if a graph G contains a (d + 1)- 
lateration graph as a spanning subgraph, then the framework (G, P) in general position also 
admits a PSD stress matrix of rank (n — d — 1). Finally, a similar result for sensor networks, 
where the graph consists of m(> d + 1) anchors is also given. 

2 The GRP and Semidefinite Programming (SDP) 

If the graph realization problem is relaxed to the problem of determining whether a realization 
of the given edge- weighted graph G exists in some unspecified Euclidean space, then this relaxed 
problem can be modeled as a semidefinite programming problem (SDP). Furthermore, one can 
find a stress matrix with the maximum rank for any given framework by solving a pair of 
semidefinite programs (see also [211 [22] HJ [10]). In particular, one can formulate a pair of dual 
SDPs where A T A is a solution to the primal problem, and the stress matrix is a solution to the 
dual problem. Here, A is the extended position matrix defined in ([1]). Next, we present one such 
formulation (for other SDP formulations of the same problem, see [Tj l5j l9[ I2T]). 

Let the inner product of two matrices R and Q be defined by R ■ Q = Trace(R T Q). An SDP 
for the relaxed graph realization problem attempts to find a symmetric matrix Y G M nxn that 
solves 

maximize • Y 

subject to (e; - ej )(e; - ej ) T ■Y = d 2 ij ,V (i< j, j) 6 E{G) (4) 
Y h 

where ej E R n is the vector of all zeros except 1 at the jth position, and Y y constrains Y to 
be symmetric PSD. A T A and P T P are both feasible solutions to Problem @ since 

(e; - ej)(ei - ejf ■ A T A = ||a* - a^ 2 = \\ Pi - Pj \\ 2 = d%, V (i < j,j) £ E{G). 

The dual of Problem ([4]) is: 

minimize £(i<,j)eE(G) w tj d^ 

subject to S := Wijifii — ej)( e i — e j) T ^ 

(i<id)eE(G) 
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Note that the dual problem is always feasible, since Wij = for all € E(G) is a feasible 

solution. In fact, this solution is also optimal, since by the weak duality theorem, is a lower 
bound on the objective. 

Prom the duality theorem, any optimal solution S of ([5]) and any feasible solution Y of 
dU) will satisfy Y ■ S = 0. This implies A T A ■ S = ASA T = 0, or AS = 0. Moreover, 
Sij = 0, V(«,j) E(G), so that any dual optimal solution is a PSD stress matrix. We say 
that the SDPs (j3J) and ([5]) admit a strictly complementary solution pair when their respective 
solutions (Y, S) satisfy rank(y) + rank(S') = n. 

Thus, the question of determining whether there is a non-trivial PSD stress matrix is equiv- 
alent to determining whether there is a non-trivial dual optimal solution, given that the primal 
problem is feasible. In particular, when the framework is universally rigid in M. d , the primal 
problem Q has a solution Y = A T A with rank (d + 1). Hence, the SDP Problems (|2J) and ([5]) 
admit a strictly complementary solution pair if and only if there is a dual optimal solution S 
for ([5]) with rank (n — d — 1). 

Proposition 1. A universally rigid framework of n vertices in W 1 , d < n — 1, always admits a 
non-trivial positive semidefinite stress matrix. 

Proof. This follows simply from Theorem 6 in [I], which states that a framework (G,P) 
in M. d admits a non-trivial PSD stress matrix if and only if there does not exist a framework 
(G,Q) in M"- 1 such that ||$ - qj\\ = \\vi~Vj\\ for all G E(G). 

□ 

The following result, stated in [21], answers the question of whether we could find a non- 
trivial PSD stress matrix, if it exists. 

Proposition 2. A primal solution Y of §J^j that has the highest possible rank among all primal 
feasible solutions, together with a dual solution S of |2U that has the highest possible rank among 
all dual optimal solutions, can be computed approximately by an SDP interior-point algorithm 
in polynomial time of n, d, and log(l/e) with error e. 

Proposition [2] also implies that if a universally rigid framework of n vertices in ~R d , d < n — 1, 
admits a rank (n — d — 1) and PSD stress matrix, then such a stress matrix can be computed 
approximately in polynomial time. However, we may not be able to compute such a stress 
matrix exactly using the SDP algorithm, even when Y is known. 

3 Main Result 

The following theorem, whose proof is given at the end of this section, is our main result. 

Theorem 3. Let (G, P) be a framework of n vertices in general position in M. d , d < n — 1, where 
G is a (d + l)-lateration graph. Then (G,P) admits a positive semidefinite stress matrix with 
rank (n — d — 1). Moreover, such a stress matrix can be computed exactly in strongly polynomial 
time, 0(n 3 + nd 3 ) arithmetic operations, if the lateration ordering and the position matrix P 
are known. 
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An n x n symmetric matrix S that satisfies condition ([2]), i.e., AS = 0, is called a pre-stress 
matrix^ Our constructive proof of Theorem [3] first generates a PSD pre-stress matrix with rank 
(n — d—1), then uses this pre-stress matrix as a basis to generate a PSD stress matrix with rank 
(n — d — 1). Recall that a stress matrix is a pre-stress matrix which also satisfies condition ([3]), 
i.e.,S ij = Q,V(i,j)?E(G). 

The following result follows from basic linear algebra. 

Proposition 3. For any framework in M. d , there exists a pre-stress matrix which is positive 
semidefinite and has rank (n — d — 1). Moreover, a universally rigid framework in M d on a 
complete graph has a rank (n — d — 1) positive semidefinite stress matrix. 

For example, the projection matrix 

i" — A T (AA T )~ l A, 

where A is the extended position matrix, is a PSD pre-stress matrix with rank (n — d — 1). 
Clearly, the projection matrix can be constructed in 0(n 3 ) arithmetic operations. 

Under the general position assumption, one can find a matrix L € ~§i nx ( n - d - 1 ) Q f the form 



L 



* : 
1 * 
1 




\0 



* * 

* * 



i * 

1/ 



that is, for k = 1, (n — d — 1), L^. 
that 



lfori = d + l + /c and = for i > (d + k + 1), such 



AL = 0, 



where A is the extended position matrix. Clearly, L has rank (n — d — 1), thus S = LL T is a PSD 
pre-stress matrix with rank (n — d—1). Such a matrix L is called a Gale matrix of framework 
(G, P) since its columns form a basis for the nullspace of A pp. 

For a (d + l)-lateration graph G with lateration ordering 1,2, ... ,n, and for a vertex k G 
V(G), let 

JV(fc) = {ie V(G) : i < k - 1 and (j, k) G £(G)}. (6) 

Thus, for such a graph, |iV(fc)| = d + 1 for each vertex k = d + 2, . . . , n. Furthermore, one can 
generate the fcth column of L, for k = 1, ...,(n — d — 1), by solving the system of linear 
equations 

Afc«i = — Od+fc+l) (7) 



1 The term pre-stress has been used by Connelly ef to mean something different, see [12J 
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where aj is ith column of the extended position matrix A, and assigning = for all i N(k). 
The above dxd linear equation system can be solved in 0(d 3 ) operations and there are n — d—1 
many of them to solve, and the formation of S takes at most 0(n 3 ) operations. Therefore, we 
have the following theorem. 

Lemma 1. The linear system ([?[) has a unique solution under the general position condition. 
Moreover, the matrix 

n—d—l 

s n = LL T = E L kLl y o 

k=l 

is a pre-stress matrix with rank {n — d — 1), and can he computed in 0(n 3 + nd 3 ) arithmetic 
operations. 

Next, we present an algorithm which uses S n of Lemma [T] as a basis to generate the desired 
stress matrix. 

3.1 A Purification Algorithm 

If the pre-stress matrix S n , as constructed in Lemma [TJ satisfies condition ([3]), i.e., Sij = 
0, V(i, j) E(G), then it is the desired stress matrix. This is true if the graph is a (d + l)-tree 
graph, that is, if there is a permutation ir of the vertices such that, 

• the first [d + 1) vertices, 7r(l), . . . , ir(d + 1), form a clique, and 

• each vertex ir(J), for j = (d + 2), . . . , n, is adjacent to the (d+ 1) vertices of a (d+ l)-clique 
in the set {vr(l), vr(2), • • • , n(J - 1)}. 

In this case, any entry in S n = LL T , for i < j and (i, j) ^ E(G), is zero. 

However, if S n is not a stress matrix, we need to zero out the entries which should be zero 
but are not, i.e., the entries Sfj ^ 0, i < j and E(G). We do this in reverse order by 
column; first, we zero out the entries Sf n ^ 0, for i < n and (i,n) ^ E(G), and then do the 
same for columns (n — 1), (n — 2), . . . , (d+3). This "purification" process will keep the pre-stress 
matrix PSD and maintain rank (n — d — 1). 

If S n is constructed from L as in the previous section, there is no need for purification of 
the last column (or row), since any entry in LL T for i <n and (i,n) E(G) is zero. But for 
general pre-stress matrices, this may not be the case. Therefore, we first show how to purify the 
last column (or row) of a PSD pre-stress matrix with rank (n — d — 1). We construct a vector 
s n <E W 1 with the elements, 

s? = -S? n ,V(i,n)<£E(G) and s n n = 1, 
and solve the following system of linear equations for the remaining entries in s n , 

i£N(n) (i,n)(£E(G) 

The right-hand-side of the equation can be formed in at most 0{nd) operations, and the dxd 
linear system can be solved in 0(d 3 ) operations. Thus, s n can be computed in at most 0(nd+d 3 ) 
operations. 
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The linear system ([8]) has a unique solution under the general position condition, and by 
construction, As n = 0. 

Lemma 2. Let S^ 1 = S n + s n (s n ) T . Then 

• AS" 1 ' 1 = 0. 

• S n ~ l >z and the rank of S n ~ l remains (n — d — 1). 

• S?- 1 = for all i < n, (i,n) E(G). 
Proof. The first statement holds, since 

AS 11 " 1 = AS n + As n {s n ) T = As n {s n ) T = 0, 

where the last step follows from the construction of s n , so that As n = 0. 

The second statement follows from = S n + s n (s n ) T h S n h 0. Thus, rank^™" 1 ) > 

rank(S" 1 ) = (n — d—1), but AS™- 1 = implies that the rank of S 71-1 is bounded above by 
(n-d-l). 

The third statement is also true by construction. In the last column (or row) of s n (s n ) T , the 
ith entry, where i ^ n and (i,n) E(G), is precisely — S™ n , i.e., 

( B nr a n\T\ _ n n _ n _ on 
\ { i )in ~ i n ~ b i ~ °in> 

so that it is canceled out in the last column (or row) of matrix S 71 ^ 1 = S n + s n (s n ) T . 

□ Note that update S" 1 " 1 = 

S n + s n (s n ) T uses 0(n 2 ) arithmetic operations. 

We continue this purification process for (n— 1), . . . ,k, . . . , (d+3). Before the kth. purification 
step, we have S k >z 0, AS k = 0, rank(5 fc ) = (n — d — 1), and 

Sg = 0, Vj > k, i < j and (i,j) E(G) 

We then construct a vector s fc G W 1 with the elements, 

»i = Si, V(i, fc) £(G), 4 = 1, and Sj fc = Va > k, 

and solve the system of linear equations for the remaining entries in s k : 

(i,k)eE(G) {i,k)£E(G) 

Again, solving this linear system takes at most 0(nd + d 3 ) operations, and by construction, we 
have As k = 0. 

Similarly, the following lemma shows results analogous to those in Lemma for the remain- 
ing columns. 

Lemma 3. Let S k ~ l = S k + s k {s k ) T . Then 

• AS 1 '- 1 = 0. 
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• S 1 y and the rank of S 1 remains (n — d — 1). 

• Sjr 1 = /or a// j>k andi< j, 

Proof. The proof of the first two statements is identical to that in Lemma [2j 
The third statement is again true by construction. Note that in the kth column (or row) of 
s k (s k ) T , the ith. entry, i > k and (i,k) E(G), is precisely —S k k , i.e., 

(s k (s k f) =s k s k = s k = -S k k , 

\ /IK 

so that it is canceled out in the kth. column (or row) of matrix S k ~ 1 = S k +s k (s k ) T . Furthermore, 
for j = (k + 1), . . . , n, the jth. column (or row) of s k (s k ) T has all zero entries, which means the 
entries in jth. column (or row) of S k ~ 1 remain unchanged from S k . 

□ 

Now we are ready to prove our main result. 
Proof of Theorem [3l 

Assume that the (d + l)-lateration graph has the lateration ordering 1,2,... , n. The matrix 
S d+2 , constructed via the process described in Lemmas [2] and El will be a PSD stress matrix 
with rank (n — d— 1), for the {d + l)-lateration graph, since after step k = (d + 3), we will have 
a matrix S d+2 that satisfies, 

AS d+2 = and S d+2 = 0, V(i, j) E(G) 

Note that the first {d + 2) vertices form a clique in G, and the principal (d + 2) x (d + 2) 
submatrix has no zero entries. This stress matrix is unique and always exists since the graph is 
a (d+ l)-lateration graph, and thus there is always a unique solution to the linear equation Q. 
Furthermore, by Lemma El S d+2 y and the rank of S d+2 remains (n — d — 1). 

There are (n — d — 2) purification steps, where each step computes a rank-one matrix s k (s k ) T 
and forms a new pre-stress matrix S k + s k (s k ) T , taking at most O(n 2 +o! 3 ) arithmetic operations. 
Thus, the computation of the max-rank PSD stress matrix uses at most 0(n 3 + nd 3 ) operations. 

□ 

We also have the following corollary: 

Corollary 1. Any universally rigid framework (G,P) in general position admits a positive 
semidefinite stress matrix with rank (n — d — 1), if G contains a (d + l)-lateration graph as a 
spanning subgraph. Moreover, such a stress matrix can be computed exactly in strongly polyno- 
mial time, C(n 3 + no! 3 ) arithmetic operations, if the lateration ordering and the position matrix 
P are known. Otherwise, such a stress matrix, together with the position matrix P, can be 
computed approximately by an SDP interior-point algorithm in time polynomial in n, d, and 
log(l/e), with error e. 

This secondary result holds because we can ignore all edges outside of the (d + l)-lateration 
spanning subgraph to prove the existence of a PSD stress matrix with rank (n — d — 1). Since 
finding a (d + l)-lateration spanning subgraph requires at least 0(n d+2 ) operations, we cannot 
actually construct such a stress matrix exactly in C(n 3 + no! 3 ) operations, if either the lateration 
ordering or the position matrix P is unknown. However, Proposition [2] implies that such a 
rank (n — d—1) and PSD stress matrix, together with the position matrix P, can be computed 
approximately in polynomial time, although not strongly polynomial. 
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4 Strong Localizability of (<i+l)-Lateration Graph with Anchors 

In this section we study the stress matrix of a sensor network, or graph localization with anchors. 
A sensor network consists of m(> d + 1) anchor points whose positions, pi,.. - ,p m G ^- d , are 
known, and n sensor points whose locations, xi,...,x n E M. d , are yet to be determined. We are 
given the Euclidean distance values d^j between p k and Xj for some (k,j), and dij between Xi 
and Xj for some i < j. Specifically, let 

N a = {(k,j) ■ dkj is specified} and N x = : i < j, dij is specified}. 

The problem is to find a realization of xi, . . . , x n E M. d such that 

\\p k - xjW 2 = & kj V(k,j)eN a 

(10) 

\\xi -xj\\ 2 = djj V{i,j)eN x . 

The semidefinite programming relaxation model for (I10p attempts to find a (d+n) x (d + n) 
symmetric matrix 

that solves the SDP 

maximize • Z 
subject to Zi : d,i:d = h 

(0; ei - ej )(0; ei - ej ) T ■ Z = d 2 {j V E N x (12) 
{-Pk]ej)(-p k ;ej) T ■ Z = (P kj V (k,j) G N a 

z y o, 



T rr I h (HI 



where (— Pfc! Cj) G M is the vector of — pk vertically concatenated with ej. Zi._^va is the d x d 
top-left principal submatrix of Z and Id is the cl-dimensional identity matrix. Z\-d^\-d = Id can 
be represented as d(d + l)/2 linear equality constraints. 
The dual of the SDP relaxation model is given by: 

minimize Id • V + + ^kjd^j 

(i,j)eN x (k,j)GN a 

subject to 5 := ( o) + 5^ ei ~ e i)(°; e * ~ e of 

+ ^ riJ} kj(-pk;ej)(-p k ;e j ) T h 0. 

Note that the dual is always feasible, since the symmetric matrix with V = 0G R dxci , variables 
Wij = for all (i, j) G N x and Wkj = for all (k,j) G iV a , is feasible for the dual. Also, each 



(13) 
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column j = (d + 1), 



,(d + n) of the dual matrix S has the structure: 




'jj 



~ J2(k,j)&N a w kjPk, 

-Wij, G N x , 

T,(i,j)EN x w ij + T,(k,j)EN a w kj- 



(14) 



Let and Wkj be called stress variables, and 5* be called a variable stress matrix for the sensor 
network localization problem. 

It is shown in [21] that both SDPs (|12p and (|13|) are feasible and solvable when there is at 
least one anchor point, the graph is connected, and there is no duality gap between the two 
SDPs. Let P be a position matrix of the n sensors satisfying constraints in (|10p . Then, the 
sensor network is said to be uniquely localizable if 



is the only matrix solution to the primal SDP ()12j) ; this is similar to the concept of universal 
rigidity. The network is said to be strongly localizable if there is an optimal dual stress matrix 
S such that 

• zs = o, 

• S + and rank(S) = n. 

It has been shown in [21 J that strong localizability implies unique localizability. 

The standard graph realization problem is equivalent to the sensor network localization 
problem without anchors; thus, the two problems are different, but closely related. For example, 
unlike in the SDP ([4]), Z constructed from A, where A is the extended position matrix, is no 
longer feasible for (|12p . although Z constructed from a position matrix P in f)15|) is feasible. 
Hence, the stresses of the dual on the anchors may not need to be balanced. As another 
example, consider a sensor network of two anchors and one sensor in M 2 , where the distances 
from the sensor to the two anchors are known. The network is not uniquely localizable, but it 
is universally rigid in graph realization, since the three points form a clique. 

However, if the sensor network has at least anchors in general position, and the graph 

realization problem has a (d + l)-point clique also in general position, then unique localizability 
is equivalent to universal rigidity, and strong localizability is equivalent to a framework on 
(n + d + 1) points having a PSD stress matrix with rank n (see |20|. 123]). The latter implies that 
the SDP pair (|12p and f)13|) admits a strictly complementary solution pair. 

Theorem 4. Take a graph G ofm(> d+1) anchor points and n sensor points with edges given 
in N x and N a , and let G be a (d + 1) -lateration graph with (P,P) in general positions. Then the 
sensor network is strongly localizable, and a rank n optimal dual stress matrix can be computed 
exactly in strongly polynomial time, 0(n 3 + nd 3 ) arithmetic operations, if the lateration ordering 
and the sensor position matrix P are known. 




(15) 



Proof. 
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We need to show that, in 0(n 3 + nd 3 ) arithmetic operations, one can compute a symmetric 
matrix S G ^(d+n)x(d+n) w j 1 j c j 1 satisfies ZS = 0, S ^ 0, rank(S') = n, and meets the structure 
condition (|14p . The proof is more complicated than that of Theorem [3j since anchor positions 
appear explicitly in the dual stress matrix. 

For simplicity and without loss of generality, we assume there are exactly (d + 1) anchors 
which are the first (d + 1) points in the lateration ordering; all other points are sensors and 
ordered 1, • • • , n. Given a position matrix P, the primal feasible solution matrix Z in (|15p can 
be written as 

,£„)-(£)«■ ") 

so that the matrix L = [-P; I n ] € M( rf + n ) xn is in the nullspace of Z, or matrix [1^ P]. 
Moreover, the matrix 

S» = LL T =( P1 *r /Ho, 
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will also be in the nullspace of Z, where rank(»S n ) = rank(L) = n. One may call S n a pre-stress 
matrix for the sensor localization problem. S n may not be a true optimal stress matrix since it 
may not meet the structure condition (I14j) . 

We now modify the elements of S n in C(n 3 ) operations, so that these conditions are main- 
tained, and the resulting matrix becomes feasible to the dual problem (|13p . Note that the 
elements in the bottom right n x n submatrix of any dual feasible S, which corresponds to the 
sensor to sensor edges, will be 



Sij 



-Wij, {hi) G N x 

,J2(i,j)£N x W ij + J2(k,j)€N a Wkj, i = j- 



Thus, for each sensor j, the jth diagonal element is the negative sum of the edge weights of w^j 
and Wkj- Since there are no constraints on V in the dual (]13p . any principal d x d submatrix 
is feasible as long as it remains positive semidefinite, and the submatrix corresponding to the 
sensor to anchor edges is feasible if, for i < d, j > d, 

Sij = - ^ W kj(Pk)i- 

(k,j)eN a 

Similar to the constructed proof of Theorem [3l while maintaining its rank n and keeping it 
PSD, we modify each column of the pre-stress matrix S n , starting with the last, (d + n), and 
continuing to column (d+1), to make it a true stress matrix, optimal for the dual. Each column 
of the matrix will be modified in the same way. 

We modify each column (d + £) of the original matrix by constructing a vector s e so that 
the the (d + £)th column of the matrix S i_1 = S e + s i (s t ) T is dual feasible. Moreover, S i ~ 1 will 
be a positive semidefinite matrix with rank n, and ZS 1 '" 1 = 0. Again, when we modify column 
d + t, the (d + j)tih column (or row) of the modified pre-stress matrix is unchanged for all j > I. 
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More precisely, when modifying the ith. column, we construct a vector s € 



pd+n 



1, sf = for i > (d + and the first (d + £ — 1) entries are 



jt ._ ( T,{k,l)eN a w k£Pk 

V Lf(i<£,e)eN x w it- e i I 



such that 



(16) 



where the (d+ 1) stress variables tD/^ and are yet to be determined, and £ K n is the vector 
of all zeros except 1 at the ith position. 

For the updated matrix S^ _1 := S £ + s^(s^) T , adding sfi(s^) T to S e will not affect any 
column (or row) to the right (or below) of column (or row) d + £. In particular, the (d + £)th 
column of S^" 1 becomes 



■e-1 

d+£ 



1 \ 

(d+e-i),(d+e) \ 



S, 



(d+i),(d+e) 



+l):{d+n),(d+£)J 



( ( -J2(k,£)GN a w kePk 



\ 



1 + S: 



(d+e),(d+e) 



\ S {d+£+l):{d+n),{d+£) ) 



In the case where £ = n, there is no last term and this becomes 



S. 



n— 1 



(( ~ E(fc,n)e7V a w knPk 
\~ Yl(i<n,n)eN x w in e i 



(d+n),(d+n) 



J 



By construction, column (d + of S i ~ 1 almost meets the the structure conditions of (JH 
To ensure S 1-1 is orthogonal to Z, or [i^ P], we determine the (d+ 1) stress variables Wk£ and 
w;^ in s e such that P]s e = 0, or equivalently, 



£-1 



5^ WHPk- Yl w uPi +P£ = S e i:d,(d+e) + /Z S d+i,(d+£)Pi- 

(k,£)eN a (i<£,£)£N x 



(17) 



1=1 



Finally to meet the diagonal entry value condition of JUJ), that is, the sum of the total edge 
stresses of sensor £ equals to the value of its diagonal element, we add 



X, w ki+ Yl Wii = 1 + S (d+e),(d+e) + Yl S d+i,(d+e)- 

(k,£)€N a (i<£,£)€N x i=£+l 



(18) 



Equations ()17p and f)18|) are exactly (d+ 1) linearly independent equations on the (d+ 1) stress 
variables, and thus there is always a unique solution. 

The modified £th column of <S^ _1 is now feasible for a dual solution stress matrix, and <S^ _1 
satisfies 



ZS l 



zs e + z*V 



o, 
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• S*- 1 = S e + s l {s l ) T hS e hO, and hence rank^ 1 ) = n. 

Repeating this process on each column, for columns [d + n), . . . , (d + 1), will result in an 
optimal dual solution matrix S° that satisfies 

• ZS° = 0, 

• S° >z 0, and rank(S' ) = n. 

That is, S° is now a true optimal dual stress matrix with rank n. Note that there are a 
total of n modification steps, and each modification step takes at most 0(n 2 + d 3 ) arithmetic 
operations. Therefore, we have proved Theorem [3} 

□ 

Similar to the secondary result for the standard graph realization problem, we have the 
following corollary. 

Corollary 2. Take a graph G of m(> d+ 1) anchor points and n sensor points with edges given 
in N x and N a , and let G contain a (d + l)-lateration spanning subgraph with (P, P) in general 
positions. Then the sensor localization problem on G is strongly localizable. Moreover, a rank 
n optimal dual stress matrix can be computed exactly in strongly polynomial time, 0(n s + ne? 3 ) 
arithmetic operations, if the lateration ordering and the position matrix P are known. Otherwise, 
such a rank n stress matrix, together with the position matrix P, can be computed approximately 
by an SDP interior-point algorithm in time polynomial in n, d, and log(l/e) ; with error e. 

The argument for Corollary [2] is analogous to that of Corollary Q3 



5 Examples 



Consider a 3-lateration framework in dimension 2 on n = 7 nodes, with ordering 1,2, ... ,7, and 
position matrix 

p _ ( -1 1 2 1 -1 -2 \ 2x7 
P ~ V 1 1 0.5 -1 -1 J GM • 

Example 1 Let 

JV(4) = {1,2,3}, N(5) = {1,3,4}, N{6) = {1, 2, 4}, N(7) = {3, 4, 5}. 
For this example, 

/ 1.5000 5.0000 -2.0000 \ 
-0.5000 3.0000 
-2.0000 -8.0000 -1.6000 
1.0000 2.0000 -2.0000 1.4000 
1.0000 -0.8000 





V o 



1.0000 

1.0000 / 
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and we have the pre-stress matrix S 7 = LL 1 



S< 



/ 31.2500 


-6.7500 


-43.0000 


15.5000 


5.0000 


-2.0000 


o \ 


-6.7500 


9.2500 


1.0000 


-6.5000 





3.0000 





-43.0000 


1.0000 


70.5600 


-20.2400 


-6.7200 





-1.6000 


15.5000 


-6.5000 


-20.2400 


10.9600 


0.8800 


-2.0000 


1.4000 


5.0000 





-6.7200 


0.8800 


1.6400 





-0.8000 


-2.0000 


3.0000 





-2.0000 





1.0000 





V o 





-1.6000 


1.4000 


-0.8000 





1.0000 / 



Note that S 7 is already a stress matrix that meets condition ([3]), so that no "purification" 
algorithm is needed. This example is not interesting, since the graph is actually a 3-tree graph. 



Example 2 Let 

AT(4) 
In this example, 



{1,2,3}, N(5) = {1,3,4}, N(6) = {2,4,5}, N(7) = {1, 3, 6}. 



L 



5.0000 


-8.0000 
2.0000 
1.0000 



and we the have the pre-stress matrix S 7 = LL T , 



S 7 



( 1.5000 
-0.5000 
-2.0000 
1.0000 



\ o 





-1.0000 


2.0000 
-2.0000 
1.0000 




-1.2500\ 


1.0000 



-0.7500 
1.0000 / 



/ 28.8125 


-0.7500 


-44.2500 


11.5000 


5.0000 


0.9375 


-1.2500X 


-0.7500 


1.2500 


1.0000 


-2.5000 


2.0000 


-1.0000 





-44.2500 


1.0000 


69.0000 


-18.0000 


-8.0000 


-0.7500 


1.0000 


11.5000 


-2.5000 


-18.0000 


9.0000 


-2.0000 


2.0000 





5.0000 


2.0000 


-8.0000 


-2.0000 


5.0000 


-2.0000 





0.9375 


-1.0000 


-0.7500 


2.0000 


-2.0000 


1.5625 


-0.7500 


V -1.2500 





1.0000 








-0.7500 


1.0000 / 



While the last column (or row) of S 7 meets condition ([3]), the rest does not satisfy (JSj) - We start 
the purification process from k = 6, where S & = S 7 . The column vector s 6 is generated by first 
assigning 

4 = -Sl 6 = -0.9375, s 6 3 = -S 3 6 6 = 0.75, s 6 6 = 1,4 = 0, 
and then solving for (s|, s|, s|) from the linear system to get 

/-0.9375\ 
-0.0625 
0.7500 
0.8750 
-1.6250 
1.0000 

V o / 
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and S 5 = S 6 + s e (s 6 ) T , 



S 5 



/ 29.6914 


-0.6914 


-44.9531 


10.6797 


6.5234 





-1.2500X 


-0.6914 


1.2539 


0.9531 


-2.5547 


2.1016 


-1.0625 





-44.9531 


0.9531 


69.5625 


-17.3438 


-9.2188 





1.0000 


10.6797 


-2.5547 


-17.3438 


9.7656 


-3.4219 


2.8750 





6.5234 


2.1016 


-9.2188 


-3.4219 


7.6406 


-3.6250 








-1.0625 





2.8750 


-3.6250 


2.5625 


-0.7500 


\ -1.2500 





1.0000 








-0.7500 


1.0000 / 



Next the column vector s 5 is generated by first assigning 



s 5 

2 



c5 
"°2,5 



-2.1016, s| = 1, s 



0, 



and solving for (sf , S3, s|) from linear system 



and S 4 = S 5 + s 5 (s 5 ) T 



( 157.4874 
-24.4489 
-230.4207 
80.8041 
17.8281 


-1.2500 



S 1 



-24.4489 
5.6705 
35.4319 

-15.5909 


-1.0625 




-230.4207 
35.4319 
338.7275 
-119.1138 
-25.6250 


1.0000 



/ 11.3047 \ 
-2.1016 
-16.4063 
6.2031 
1.0000 





80.8041 
-15.5909 
-119.1138 
48.2444 
2.7813 
2.8750 




17.8281 


-25.6250 
2.7813 
8.6406 

-3.6250 




One can see that S is now a desired stress matrix for Example 2. 





-1.0625 


2.8750 
-3.6250 

2.5625 
-0.7500 



-1.2500\ 


1.0000 



-0.7500 
1.0000 / 
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